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We developed the deformed quasi-particle random phase approximation (DQRPA) 
to describe various properties of deformed nuclei and applied to the evaluation of the 
Gamow-Teller (GT) transition strength distributions which data can be extracted 
from the charge exchange reactions (CEXR). Our calculations started with the single- 
particle states calculated by the deformed axially symmetric Woods-Saxon poten- 
tial. Pairing correlations of nucleons, neutron-neutron and proton-proton as well 
as neutron-proton pairing correlations, are explicitly taken into account at the de- 
formed Bardeen Cooper Schriffer (BCS) theory leading to the quasi-particle concept. 
The ground state correlations and the two-particle and two-hole mixing states are 
included in the deformed QRPA. In this work, we use a realistic two-body interaction 
given by the Brueckner G-matrix based on the Bonn potential to reduce the ambi- 
guity on the nucleon-nucleon interactions inside nuclei. We applied our formalism to 
the GT transition strengths for 76 Ge, 82 Se, 90 Zr, and 92 Zr, and compared to available 
experimental data. The GT strength distributions turn out to be sensitive on the 
deformation parameter. We suggest most probable deformation parameters for the 
nuclei by adjusting GT strength distributions to the experimental data and the Ikeda 
sum rule, which sum rule is usually thought to be satisfied more or less even under 
the one-body current approximation without introducing the quenching factor, if 
high-lying excited states are properly taken into account as in our approach. Most 
of the experimental GT strength data of the nuclei, in particular, on the high-lying 
excited states beyond one nucleon threshold are confirmed in our DQRPA by the 
deformation parameter. 
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I. INTRODUCTION 

In the core collapsing supernovae (SNe), medium and heavy elements are believed to be 
produced by the rapid and slow successive neutron capture reaction, dubbed as r-process and 
s-process, respectively. In these processes, many unstable neutron-rich nuclei are produced 
iteratively and decay to more stable nuclei at their turning points in the nuclear chart. These 
r- and s-processes play vital roles of understanding the nuclear abundances in the cosmos 

Since most of the nuclei produced in the processes are thought to be more or less deformed, 
we need to explicitly consider effects by the deformation in the nuclear structure and their 
reactions in the network calculations of the processes. One interesting process related to the 
deformed nuclei may be the rapid proton process (p-process), which is thought to be occurred 
on the binary star system composed of a neutron star and a companion star, usually white 
dwarf. Because of the strong gravitation on the neutron star surface, one expects hydrogen 
rich mass flow from the companion star. Since the high density and low temperature on 
the neutron star crust makes electrons degenerated, beta decays are usually Pauli blocked 
by the degeneration, which gives rise to the pycno- reactions [2]. It means that neutron rich 
nuclei may become stable nuclei, while unstable nuclei may become stable in this situation. 

Up to now, main conventional approach to understand the nuclear structure is based 
on the spherical symmetry. In order to describe neutron-rich nuclei and their relevant nu- 
clear reactions occurred in the nuclearprocesses, one needs to develop theoretical formalism 
including explicitly the deformation Ref. exploited the Nilsson basis to the de- 

formed quasi-particle random phase approximation (DQRPA). But the two-body interaction 
was derived from the effective separable force. The realistic two-body interaction derived 

n 

from the realistic force is firstly exploited at Ref. |4| only with neutron-neutron (nn) and 
proton-proton (pp) pairing interactions which have only isospin T = 1 and J = interaction. 
But, to properly describe the deformed nuclei, the T= and J = 1 pairing should be also 
taken into account because it is intimately related to the deformation by the J = 1 pairing. 

In this work, we extend our previous QRPA based on the spherical symmetry 5| to the 
DQRPA. The spherical QRPA has been exploited as a useful framework for describing the 
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nuclear reactions sensitive on the nuclear structure of medium-heavy and heavy nuclei 
where, as the mass number increases, the application of the shell model may have actual 
limits because of tremendous increase of configuration mixing due to the deformation. This 
paper is organized as follows. In Sec. II, we introduce detailed formalism for the DQRPA and 
the Gamow- Teller (GT) strength. Numerical results and related discussions are presented 
at Sec. III. Summary and conclusions are addressed at Sec. IV. 



II. FORMALISM 
A. Total Hamiltonian 



We start from the following nuclear Hamiltonian 

H = H + H int , (1) 

H o= e ap a a>c< apaQ ,c apcia i (a =p,n) , 

ap a a' 

Hint = Vapaa'f3p l3 l3'fp~ / 'y'6p s 6'cl ipaa ,C^ pi3!3/ CSp s siC 1 p i y, 
a/3-ySpapppyps, a'fi'-y'S' 

where the interaction matrix V is the anti-symmetrized interaction with the 
Baranger Hamiltonian [7] in which two ^ factors, from J and T coupling, are 
included, so that the H int in Eq. (TTJ is equivalent to the usual H mt = 

lEa^p^p,, a'i3' 1 'S'^ r c l p a a'i3p^' 7 p^'5p s 6'Cap ol a' c ^ppP' c Sp s S'C'yp^'- Th e Greek letters denote 
proton or neutron single particle states with the projection Q of the total angular mo- 
mentum on the nuclear symmetry axis and the parity n. The projection Q is treated as 
the only good quantum number in the deformed basis. The p a (p a = ±1) is the sign of 
the angular momentum projection Q of the a state. The isospin of particles is denoted as 
a Greek letter with prime, while the isospin of quasi-particles is expressed as a Greek letter 
with double prime as shown later. 

Therefore, the operator c' a , (c apaa ') in Eq.(l) stands for the usual creation (destruction) 
operator of the real particle in the state of ap a with the angular momentum projection Q a 
and an isospin a', and c Bipa a'i= ( — l) Ja+ma, c_ apaQ ,/) is the time reversed operator of c apaCt '. 
Namely our intrinsic states are twofold-degenerate, i.e. fl a state and its time-reversed state 
— Q a . e apaQ i means the single particle state energies. 
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In cylindrical coordinates, eigenfunctions of a s.p. state and its time-reversed state in 
deformed Woods-Saxon potential are expressed as follows 

\a Pa = +1 >= Yl^NnA \ N > °<* = A « + V2 > (2) 

Nn z 

+ b( hlnJ N , n z ,A a + l,Q a = A a + 1-1/2 >}, 
\ap a = -1 >= \ N > Uz > _Aa ' Qa = ~ Aa ~ l l 2 > 

Nn z 

-b ( N ln a n z , -A a -l,Q a = -A a -l + 1/2 >], 

where N = n±_ + n z (n± = 2n p + A ) is a major shell number, and n z and n p are numbers of 
nodes on the deformed harmonic oscillator wave functions in z and p direction, respectively. 
A is the projection of the orbital angular momentum onto the nuclear symmetric axis z. The 
coefficients 6^2 n anc ^ "nIq are obtained by the eigenvalue equation of the total Hamiltonian 
in the Nilsson basis. The 2nd terms in Eg. (j5J) have the same projection Q a value as the 
1st term, but retain another orbital angular momentum because of a flipped spin. Particle 
model spaces we used depend on the nuclei as shown later on. Single particle spectrum 
obtained by the deformed Woods Saxon potential is sensitive on the deformation parameter 
^2 defined as 

R{9) = R,{l + f5 2 Y w {9)) , (3) 

where Rq = 1.2A 1//3 fm for the sharp-cut radius Rq [8| and Y 2 o is spherical harmonics. The 
customary parameter e = 3(u± — u 3 )/(2uj± + w 3 ) used in the deformed harmonic oscillator 
is related to as 0% ~ (2/3)a/47t/5 e at the leading order. We transform the Hamiltonian 
represented by real particles in Eq. (TjQ) to the quasi-particle representation through the 
following Hartree Fock Bogoliubob (HFB) transformation, 

The Hamiltonian, then, can be expressed in terms of the quasi-particles as follows 

H = H Q + E aa na i,a aa n + H qp , in t . (5) 

'act" 

Using the transformation of Eq. (jl]), finally we obtain the following HFB equation 
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where E aa " is the energy of a quasi-particle with the isospin quantum number a" in the 
state a. If we neglect A np , this equation reduces to the standard Deformed BCS (DBCS) 
equation. 



B. Spherical and deformed wave functions for a single particle state 

Since various mathematical theorems regarding quantum numbers may not be easily used 
in the deformed basis, it is more convenient to play in the spherical basis. In addition, we 
exploit the G-matrix based on the Bonn potential in order to pin down plausible ambiguities 
on the nucleon-nucleon interaction inside a deformed nucleus. Since the G-matrix is calcu- 
lated on the spherical basis, we need to represent the G-matrix in terms of the deformed 
basis. Here we present as to how to transform deformed wave functions to spherical wave 
functions. 

The deformed harmonic oscillator wave function, \Nn z A a Q a (= A a + E) >= liVn^Ao, > 
|S > in Eq. (T5]) can be expanded in terms of the spherical harmonic oscillator wave function 
\NqIA > |S > as follows 

AT I — N Q~ l 

\Nn z A a > |E>= E A NnJ r ~~ \NvlA a > |E >, (7) 

N =N,N±2,N±4,... i=JV ,iVo-2,JVo-4,... 

|iV /A a > |E >= ^C^\jN lj Q a > . 

j 

Here the spatial overlap integral A™ a =< NolA\Nn z A > is calculated numerically in the 
spherical coordinate system. Roman letter indicates quantum numbers (n r lj) of a nucleon 
state, where n r is the radial quantum number and / is the orbital angular momentum. C ; J A a is 
is the Clebsch-Gordan coefficient of the coupling of the orbital and spin angular momenta 
to total angular momentum with the projection Q a . Therefore, the expansion can be simply 
written as 

\an a >=Y,B: \an a >, (8) 

o 
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with the expansion coefficient summarized as B% = ^2 Nnz % C^i^jvn^A ^Nn z T,- 

The DBCS Eq. fl6]) is solved by using the Brueckner G-matrix calculated from the realistic 
Bonn CD potential for the nucleon-nucleon interaction in the following way. The pairing 
potentials A p , A n and A pn in Eq. ([6]) are calculated as 

& ap « P = ~\ (2q + X 1)1/2 9^M a n)FiaaaF'ncG{aacc, J) (2c + l) 1/2 {u* lp v lpc + u* p v 2pc ) , 

J,c 

(9) 

where = 5" £>f C!j^n a j^n {K' = Q a + fya) is introduced for the transformation to the 

deformed basis of G-matrix. Here K', which is a projection number of the total angular 
momentum J onto the z axis, is selected K' = at the BCS stage because we consider 
the pairings of the quasi-particles at a and a states. A anSn is the same as Eq. (jUJ) with 
replacement of n by p. In order to renormalize the G-matrix, strength parameters, (7p air , 
(?p air and are multiplied to the G-matrix [5| by adjusting the pairing potentials to the 
empirical pairing potentials, A p mp , A^ mp and A p ™ p . The empirical pairing potentials of 
protons (neutrons) and neutron-proton are evaluated by the following symmetric five term 
formula for the neighboring nuclei 

A p cmp = -[M(Z + 2,N)- AM(Z + 1, N) + 6M(Z, N) (10) 
8 

-AM(Z - 1, N) + M(Z - 2, N)\ , 

A n omp = -[M(Z, N + 2) - 4M(Z, N + 1) + 6M(Z, N) (11) 
8 

-4M(Z, N-l) + M(Z, N-2)], 

A pn cmp = ±j{2[M(Z,N + l) + M(Z,N-l) + M(Z -1,N)+M(Z + 1,N)} (12) 
-[M(Z + 1, N + 1) + M(Z - 1, JV + 1) + Af (Z - 1, JV - 1) + M(Z + 1,N— 1)] 
-4M(Z,iV)} , 

where +(-) stands for even(odd) mass nuclei. As for masses in Eqs. (10) ~ (12), we use 
empirical masses. For the neutron-rich nuclei far from the stability region, the masses are 
still unknown. We, thus, are going to use the values from the extrapolation or the values 
derived from the liquid-drop model, A w 12A _1//2 MeV. More exact mass measurement of the 
exotic nuclei from the penning trap [h| are greatly helpful for these empirical parameters. 
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C. Description of an excited state by the DQRPA 



We take the ground state of a target even-even nucleus as the DBCS vacuum for a quasi- 
particle. In the following, we show how to generate an excited state in a deformed compound 
nucleus. Since the deformed nuclei have two different frames, the laboratory and the intrinsic 
frame, we need to consider the relationship of the two frames. The GT excited state in the 
intrinsic frame of even-even nuclei, which is described by operating a phonon operator to the 
BCS vacuum Q m K \QRPA >, can be transformed to the wave function in the laboratory 
frame by using the Wigner function 0,ip) as follows 

\lM(K),m> = ^[V 1 M _ K (<j ) ,e^)Ql,K\QRPA> (for K — 0), (13) 



\lM(K),m> = ^—- 2 [V l MK {^9^)Ql,K 

+ (-l) 1+K V 1 M „ K ( ( j ) ,6^)Ql_ K \QRPA > (for K = ±1), 

where M(K) is a projection of the total angular momentum onto the z (the nuclear symme- 
try) axis. Of course, the K is accepted as only a good quantum number in deformed nuclei. 
The QRPA phonon creation operator Qp m K acting on the ground state is given as 

Ql,K = E [ X ^"W>K)AKaa"(3{3"K) - ^ K) A(aa" 'ffl> 'K)\. (14) 

aa"/3l3" 

with pairing creation and annihilation operators comprising two quasi-particles defined as 

A^(aa"(3f3"K) = a ] aall a^ in A(aa"(3f3"K) = a^a^,,, (15) 

where bar denotes a time-reversal state for a given state. The quasi-particle pairs in two- 
particle states, a and j3, are chosen by the selection rules Q a — = K and ir a irp = 1. 

Two-body wave functions in deformed basis are calculated from the spherical basis as 
follows 

\aP>=J2FL%\ab,JK> : (16) 

abJ 

where 2-body wave function in the spherical basis is given as \ab, JK >= 
^2jCjj^ a j a \an a > \bfl b >, and the transformation coefficient is calculated as F^ b = 
Bl B%(-l) j e- n P C^ aj a which has a phase factor (-l)w-ty arising from the time- 
reversed states ft. B% is defined below Eq.(8). 
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Our QRPA equation in deformed basis leads to the following form. 
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where 1 and 2 denote the isospin quantum numbers of protons and neutrons. The amplitudes 
X aa "pp" an d Y a a"pp", which stand for forward and backward going amplitudes from state 
aa" to f3/3", are related to X a a"pi3" — v^&aa" pp" X aa " pp" and _Y a a"pp" — y^^aa" pp"Yaa" pp" ■ 
Caa" pp" — 1 if Oi — ft and a" = (3", otherwise a aa "pp" = \/2 [5j. The A and B matrices are 
given by 



(E aa ii + -E'/3^'')5a7^a"7"^/35^/3"<5" ~ ° aa" PP" a ^" 55" 
x \9vv^ U aa" u pp" u -fi' u 55" + v aa"Vp-^'V y yiVg S „) V^ ^g 
+ 9phi, v >aa" v pp>' v ' r i r i" v 55" + v aa" u pp" v 'yy"' lJ 'Sgr') ^ 

+ 9 P h(u 



(18) 



7" 5" 

o,9, 7<5 



(if) = — 0~aa"PP" Cr - r y"55" 

x [ — 9pp{Uaa"UppiiV^nVggii + V a a"V p^iiUyy/'Uggn) g 

+ gph(Uaa"Vpp"V~ r y"Ussr> + V aa "Upp"U'ry"V56") K«5, 7/3 

+ 9ph{ u aa" v pp" U 'Y'r" v 5S" + v aa"Upp"V~ry"UsS") Vary, 5p]i 



(19) 
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V 5, js = J2J2 F ^ F i^G(abcd, J) , (20) 

J abed 
J abed 

v £,6fi = J2J2 F ^ c F^ G (acdb, J). 



J abed 

where F JK ^, = B® B? C/'o' a n with non-zero K' = Q a + Qg, and u and v coefficients are 
determined from deformed HFB calculation with the pairing strength g™ aiT , and g^ ir 
adjusted to the empirical pairing gaps A nn , A pp and A np , respectively. E aa » indicates the 
quasi-particle energy of the state a with the quasi-particle isospin a". The G(F) matrices 
are two body particle - particle (hole) matrix elements obtained as solutions of the Bethe 



- Goldstone equation from the Bonn potential 



□ 



If we do not consider the np pairing 



correlations, the A and B matrices can be expressed in the following simple form, which is 



the same as the those of Saleh 

^3,75 = ( E <xp + Epn) S av$pp'$p5 S nn> (21) 

~~ 2 [g pp (u a pUp n u rp iu$ n , + v ap Vp n Vj P /Vg n ,) V a ^ ^ 

+ 9ph( u ap v pn u JP ,v Sn' + v apUj3 n V~ /p /U$ n ,) V aS ^ 7(9 ], 
B P a$, 7<F (-^0 = ~ ^ [~9pp( u apU(j n V lp iV$ n , + V ap V p n U lp iU$ n ,) ^ (22) 

+ 9ph{u a pVp n v lp iu$ n , + v ap Up n u lp iv^ n ,) V aSi 7(9 ] . 

D. Description of Gamow- Teller Transition 

The /3 ± decay operator, /9 X , is defined in the intrinsic frame as 

fcn = Yl < a pPc\ T+ °K\PnPp > cjc n , /3+ = 0ij = {~YK-^ (23) 

OtppafinPp 

in which the pf 1 transition operators are related with those in the laboratory system 
operator as follows 

& = E D ^Wi- (24) 
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The ^ transition amplitudes from the ground state of an initial nucleus to the excited state, 
the one phonon state in a final nucleus, are expressed by 

< l(K),m\Px\ QRPA > (25) 

= ■N" aa "p a pp"p p < Old 'pp a \o '(3' ' Tip p > [u paa "V n ppri X aa ii ppn t K + V paa iiU n ppiiY aa ii ppu j<] 

< l(K),m\0£\ QRPA > 

= Naa"f3p" < Ota"pp a \<T K \[3(3"npp > [u paa "V n pp"Y aa " pp" >K + V paa "U n pp"X aa 'i /3i3",k} 

aa"p a /3p"p l3 

where | QRPA > denotes the correlated QRPA ground state in the intrinsic frame and the 
nomalization factor is given as N aa "pp"{J) = a/1 — S a pS a "p"(— 1) J+T /(l + 8 a p8 a "p")- The 
Wigner functions are disappeared by using the orthogonality of two Wigner functions from 
the operator and the excited state, respectively. This form is also easily reduced to the 
results by proton-neutron DQRPA without the np pairing 

< 1 (K),m \p-\ QRPA> (26) 

= Yl < a PPa\ T+a K\PnPl3 > [Ua p Vp n X ap p n:K + V ap Up n Y ap p ntK ] , 

OLppaPnPp 

< l(K),m\fy\ QRPA > 

= ^ < a p p a \T + a K \/3 n pp > [u ap Vp n Y ap p n>K + v ap up n X ap p niK ] . 

OLppaftnPp 

Tere single particle matrix elements of {a p p a \T + a x\Pnp p) can be expressed in deformed basis 
< a p p a \T + a K=0 \(3 n pp >= 5 np n nPa ^[^fl p ^ A - fefeLn^fc^nJ: ( 27 ) 

Nn z 

< a P p a \r + (T K=1 \(3 n pp > = -V25 npnn+1 J^^nL^nL^ = PP = + 1 ) ( 28 ) 

Nn z 

= +V25n p n n +i ^ b^l^b^l^ (p a = Pp = -1) 

Nn z 

= -V2^i^ n -i J2 b Nln p b Nln n (Pa = +1, Pp = ~l), 

Nn z 

< a p p a \r + a K= . 1 \f3 n pp > = V25n p n n -i ^ ^in p & !vin n (p* = Pp = + 1 ) ( 29 ) 

Nn z 

= -V2S np a n -i fi,n P 6 fel,n„ (p« = Pp= ~ l ) 

Nn z 



2 °° n 2 

Nn z 
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To compare our theoretical results to the experimental data, the GT(=p) strength functions 
and their running sums (total strengths) are calculated as 



B GT {m) = I < iW'HI&ll QRPA > | 2 , (30) 

K=0,±1 

P>gt(™) = \<1{K)M\Pk\\QRPA>\ 2 , 

K=0,±1 

S gt= E ^|<l(ir),m||/3-||g J RPA>| 2 , (31) 

K=0,±1 m 

S GT= E J2\ <1 ( K ^ m ^K\\QRPA>\ 2 . 
K=0,±1 m 

E. Ikeda Sum Rule 

Numerical results for total GT(±) strengths, Sq T and Sq T , in Eq. fl3Tj) are investigated 
through the Ikeda sum-rule (ISR), which is known to be satisfied more or less independently 
of the constructed excited states by any nuclear models, 

S GT - S GT = 3( N ~ Z )- ( 32 ) 

The ISR within the Deformed QRPA, denoted as ISR II, is given by 

(S G rp — Sq T )j SR jj 

= E Et I < 1( ^ K )M I QRPA > | 2 - I <l(K),m\ 1 1 QRPA> | 2 ] 
isr=o,±i m 

a p pa/3nPj3 K=0,±1 m 

(33) 

If we use the closure relation for the excited states, the ISR which we denote as ISR I is 
shown to be easily calculated as follows 

(Sgt ~ Sgt)isr i 

= E E \<ot p n p \T + a K \P n n n >\ 2 (vl-v 2 ap ) . (34) 

CtpPaPnPp K=0,±1 

Since we use the limited particle model space in deformed basis under the one-body current 
without the A excitation, the above sum rule is a bit broken, but it may be used to test our 
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nuclear model and their numerical calculations as shown in our numerical results. On the 
other hand, the ISR in spherical basis is easily shown to satisfy the ISR as follows 



S gt~ S gt =$^l <a p \T + a\b n > \ 2 {v 2 n -v 2 p ) (35) 



2 

1 1 1 " 



J2(2ja + l)(2j b + l)5 nanb 5 lalb l 2 2 \ {vl-v 2 p ) 

a p b„ I jb ja la J 

3 ^(2j b + l)v 2 n - 3 5^(2 Ja + l)*; 2 = 3(N - Z). 



III. RESULTS 

We calculated the Gamow-Teller (GT) strength, B% T in Eq. (30), for 76 Ge, 82 Se, 90 Zr, 
and 92 Zr within the DQRPA. Those nuclei are selected to represent medium-heavy nuclei 
which have experimental data relevant to the GT strength. The single particle states are 
used up to 4hu for all nuclei, in the spherical limit. 

Since the GT strength distributions turn out to be sensitive on the deformation parameter, 
/?2 in Eq.(3), we exploited various values of the deformation parameter, j/^l < 0.3, as the 
default values. For the pairing interaction, the strength parameters g™^ and g^ &h in Eq.(9), 
which are introduced to renormalize the finite Hilbert model particle space, are adjusted to 
reproduce the empirical pairing gaps through symmetric five term formula Eqs.(10)-(12) jsj]. 
Corresponding values of the deformation parameter 02, theoretical and empirical pairing 
gaps are tabulated in Table I with the ISR I and II in Eqs.(33) and (34). It is a remarkable 
point that theoretical ISR results are well satisfied within 15 % maximally, and nearly 
independent of the deformation parameter, /3 2 , as required in the sum rule. 

The single particle state energies adopted from the deformed Woods Saxon potential nat- 
urally depends on the parameter /3 2 . The deformation of nuclei may be conjectured to come 
from macroscopic phenomena, for example, the core polarization, the high spin states and 
so on. Microscopic reasons may be traced to the tensor force in nucleon-nucleon interaction, 
which is known to account for the shell evolution according to the recent systematic shell 
model calculations 12|, [l3j]. For example, T = 0, J = 1 pairing, which is associated with 
the 3 S"o tensor force, may lead to the deformation contrary to the spherical T = 1, J = 
pairing. Therefore, the deformation parameters adopted in this work may include implicitly 
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and effectively such effects, because the single particle states from the deformed Wood Saxon 



potential show a strong dependence on the fa |l4j . 

The deformation parameter fa may help us to conjecture the nuclear shape through the 
intrinsic quadrupole moment, Q = a/ 167r /5(3 / Air) ARq fa, where Rq = 1.2y4 1//3 fm for the 
sharp-cut radius Rq [15I, Q|- In experimental side, the deformation parameter fa can be 
extracted from the E2 transition probability, Q = ^/l67iB(E2)/5e 2 . Since we do not have 
yet enough E2 data to be exploited in the nuclei considered here to our knowledge, we refer 
to theoretical results by the relativistic mean field theory (RMF) 17J. But the sign of the 
fa is sometimes still uncertain and also the coexistence of prolate and oblate shapes are also 
possible. In this work, such coexistence and fa deformation are not taken into account. 

In the DQRPA stage, we took the particle-hole and particle-particle strength parameters, 
g p h and g pp in Eqs.(18)-(19) as 1.0 for for the four nuclei. Actually g p h might be determined 
from the Gamow Teller Giant Resonance (GTGR), while fine tuning of g pp are usually 
performed for the double beta decay j^]. In this work, we fix them because the GT strength 
distributions are not so sensitive on the values. 

Since the nuclei considered here are expected to have large energy gaps between proton 
and neutron spaces, we considered only the nn and pp pairing correlation although the 
formalism is presented generally. For example, in the neutron-rich nuclei of importance in 
the r-process, the np pairing may not contribute so much. But for the p-process, the np 
pairing could be more important than in the neutron-rich nuclei because of the adjacent 
energy gaps of protons and neutrons. The calculations for the neutron-deficient nuclei in 
p-process are in progress by the explicit inclusion of the np pairing correlations. 

In the following, we show the GT strength distribution in terms of the fa parameter, as 
a function of excited energy of parent nucleus. Therefore, experimental data, which were 
usually measured w.r.t. the the ground state of daughter nuclei, are presented by subtracting 
the empirical Q values from the measured data. 

In Figsfj]and[2l we show the GT(-) strength of 76 Ge as a function of the excitation energy 
E ex w.r.t. the ground state of 76 Ge, whose Q value is 0.9233 MeV. The uppermost two panels 
are the experimental data from the 76 Ge(p,n) 76 As reaction at 134.4 MeV 18], which show 
a strong GT state peak around 12 MeV. FigJTJ represent the calculated GT strength with 
prolate shapes, fa = 0.1 ~ 0.3, and FigfSJis for oblate shapes, fa = —0.1 ~ —0.35. The GT 
strength distributions are widely scattered owing to deformation. 
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In particular, the results for /3 2 = 0.35, in which the ISR I and II are almost satisfied as 
shown in Table 1, nicely reproduce the peak on the high- lying excited GT states around 12 
MeV. Therefore, the redistribution of the GT strength by the deformation is not a simple 
wide scattering, but should be understood as the redistribution well matched to the high- 
lying excited states to be confirmed in the experiments. Our /3 2 value is consistent with 
the prolate shape suggested by the RMF calculation, whose value was /3 2 = 0.157 [r3] a bit 
smaller than the present value. 

The GT(-) strength, B(GT— ), on 82 Se is presented, in Figj3j as a function of the ex- 
citation energy E ex w.r.t. the ground state of 82 Se (whose Q value is 0.092 MeV) for 
02 — ±0.1 ~ ±0.2. The uppermost panel of Figj3] is the experimental data from the 
82 Se(p,n) 82 Br reaction at 134.4 MeV 18j. Particularly, the strong strength on the high-lying 
excited GT state around 12 MeV is neatly reproduced by the /3 2 = 0.2 at (c) panel. For a 
reference, the /3 2 value from the RMF is 0.133 Our theoretical calculations address a 
possibility of another peak around 16 MeV. 

The forward GT amplitudes X^/ppn^ for three highest GT transition peaks at 10.74 
and 11.51 MeV in (c) panel of Figj3] are shown in FigJH Main transitions turn out to 
come from the (431 3/2) (440 1/2) at 10.74 MeV, and (301 3/2) (301 1/2) and 
(411 3/2) — > (411 1/2) transition at 11.51 MeV around the Fermi surface, respectively. They 
are thought to result from the 2p-2h configuration mixing stemming from the deformation 
around the Fermi surface. 

Figsj5] and [6] show the GT(=p) strength on 90 Zr as a function of the excitation energy 
E ex w.r.t. 90 Zr, whose Q values for GT(=p) are 6.111 and 2.280 MeV, respectively. The 
uppermost panels (a) stand for the experimental data on GT(-) and GT(+) deduced from 
the 90 Zr(p,n) 90 Nb reaction 90 Zr(n,p) 90 Y at 293 MeV fl9|. respectively. Panels (c)~(d) are 
the results by the DQRPA for two different prolate and oblate 02 values. For a reference, 
in panel (b), we show results by the spherical QRPA {5, Q]. The ISR I and II are almost 
satisfied at each /3 2 value for both B(GT(=p)). If we recollect that the experimental GT 
strengths on the high-lying states beyond 30 MeV actually include the contributions by the 
isovector spin monopole (IVSM) [3] which are not considered in the present calculation, our 
spherical QRPA (SQRPA) (b) and DQRPA (c) for /3 2 = 0.1 results seem to be consistent 
with the data. 

It means that the 90 Zr is thought to be near to a spherical shape, i.e. almost spherical. 
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But the DQRPA results by (3 2 = 0.1 show peaks below Qec — 6.111 MeV of 90 Zr which are 
inconsistent with the data. Since our DQRPA is based on the intermediate deformation, 
the extension to small 02 values may not be a proper treatment, if we understand that the 
angular momentum projection Qj may have angular momenta higher than j, and may go 
back to different j values in the /?2 = limit. One more interesting point is that almost 
spherical nuclei, such as 90 Zr, did not show the high-lying GT(-) excited states appeared 
on 76 Ge and 82 Se. The high-lying states at the GT(+) transition is also neatly explained as 
shown at (c) in FigJBl 

FigJTJ show the GT(-) strength of 92 Zr as a function of the excitation energy E ex w.r.t. 
the ground state of 92 Zr, whose Q value is 2.005 MeV. The uppermost panel (a) of FigJT] 
represent the experimental B(GT~) values extracted from 92 Zr(p,n) 92 Nb reaction at 26MeV, 
which energy was too low to expect high-lying excited states. Excitation energies observed 
in this reaction are located below 9 MeV. Panels (b) is the our results for (3 2 = 0.1. 

IV. SUMMARY AND CONCLUSION 

To describe deformed nuclei, we performed the deformed axially symmetric Woods-Saxon 
potential, the deformed BCS, and the deformed QRPA with realistic two-body interaction 
calculated by Brueckner G-matrix based on Bonn CD potential. Results of the Gamow- Teller 
strength, B(GT±), for 76 Ge, 82 Se, 90 Zr, and 92 Zr show that the deformation effect leads to a 
fragmentation of the GT strength into high- lying GT excited states, in particular, on 76 Ge 
and 82 Se which GT data were already measured at the charge exchange reaction experiments 
by the 134.4 MeV proton beam. Those states are shown to be properly explained by the 
deformation effects. 

These high-lying excited GT states may affect seriously relevant nuclear reactions, partic- 
ularly for neutrino-induced reactions, exploited in the nucleosynthesis because the emitted 
neutrinos from the proto-neutron star may have a high energy tail up to tens of MeV energy 
range. Since the experiments relate to the neutrino reaction on such a high energy would be 
a very challenging task in the present neutrino factories, the extraction of the high-lying GT 
states from charge exchange reactions could be very useful for understanding the neutrino 
reaction in the cosmos, if we recollect that the GT transitions are main components for the 
neutrino-induced reaction. 
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More systematic analysis of deformed nuclei by our DQRPA are under progress by in- 
cluding light nuclei characterized by the shell evolution or the inversion island. In the light 
nuclei or neutron-deficient nuclei, which may have small energy gaps between protons and 
neutrons and small N-Z values, the neutron-proton pairing correlations could definitely af- 
fect nuclear (5 + decays and Ikeda sum rule. Therefore GT(+) strength distribution data on 
the nuclei could be a stringent test of nuclear models. 
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TABLE I: Deformation parameters and empirical (theoretical) pairing gap parameters Acm 
(A^ 11 ) used in this work. The ISR in the last column denotes the Ikeda sum rules I (Eq.(34)) 
and II (Eq.(33)) as a percentage ratio to 3(iV — Z). The particle-particle (particle-hole) strength 
parameters are exploited as g pp = 1.0 (g p h = 1.0) for all nuclei. 

nucleus /3 2 AP m (MeV) A t p h (MeV) A£ m (MeV) A? h (MeV) ISR I, II(%) 



0.1 1.575 1.538 97.04, 106.03 

0.2 1.578 1.540 96.94, 102.40 



76p p 


yj.oo 


1 ^RzL 


1 ^fi9 


1 

l.OOO 


1 ^A^ 






-0.1 




1.578 




1.542 


97.20, 102.74 




-0.2 




1.561 




1.549 


97.04, 109.01 




-0.3 




1.566 




1.562 


97.34, 115.27 




0.1 




1.502 




1.579 


96.51, 106.45 


82 Se 


0.2 


1.409 


1.419 


1.544 


1.547 


95.62, 107.55 




-0.1 




1.410 




1.564 


96.58, 104.89 




-0.2 




1.568 




1.549 


96.57, 109.69 


90 Zr 


0.1 


1.247 


1.267 


1.705 


1.713 


97.38, 107.41 




-0.1 




1.269 




1.714 


97.88, 88.40 


92 Zr 


0.1 


1.357 


1.378 


0.841 


0.847 


97.32, 98.18 
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FIG. 1: (Color online) Gamow- Teller strength distributions B(GT~) on 76 Ge as a function of the 
excitation energy E ex w.r.t. the ground state of 76 Ge. Experimental data denoted as filled (red) 
points in the uppermost panels are deduced from the 76 Ge(p,n) reaction at 134.4 MeV [1.8]. In 
each panel, we indicate @2 = 0.1 ~ 0.35 for prolate shapes. 
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2: (Color online) The same as in Fig(TJ but for oblate shapes, $2 = —0.1 ~ — 0. 
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FIG. 3: (Color online) Gamow- Teller strength distributions B(GT ) on 82 Se as a function of the 
excitation energy E ex w.r.t. the ground state of 82 Se. Experimental data denoted as filled (red) 
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FIG. 4: Forward amplitudes X of the two highest GT strengths in (c) panel of Figj3] at 10.74 
and 11.51 MeV w.r.t. the ground state of 82 Se among the GT transition configuration mixing. 
Transition (4 3 1 3/2, 4 4 1/2) is the only physical component to the GT state at 10.74 MeV, 
and (3 1 3/2, 3 1 1/2) and (4 11 3/2, 4 11 1/2) transitions are main components for the 
GT state at 11.51 MeV. 



22 



4h 90 Zr(p J n) 90 Nb l 



CD 2 |- 
m 





(a) Exp - 



10 







O 5 - 

CD 









— 4 - 

CD 









4 - 

h- 

O 2 - 
CD 





,*, i*. ... i ,',«,t,«i», V.Vi',',*,',' 



10 20 30 40 50 
i 1 1 



Mi 



T 



T 



T 



T 



T 



(b) QRPA - 



10 20 30 40 50 

r-j— i— i— r 



(Si 



4%i . ■ ■ ■ 



I I I I 



10 20 30 40 50 

r-j— i— i— r 



(d) 7= - 0.1 







i i i | i i i i | r- 

10 20 30 
E ex [MeV] 



40 



50 



FIG. 5: (Color online) Gamow- Teller strength distributions B(GT~) on 90 Zr as a function of the 
excitation energy E ex w.r.t. the ground state of 90 Zr. Experimental data denoted as filled (red) 
points in the uppermost panels are deduced from the 90 Zr(p,n) reaction at 293 MeV [19j]. In each 
panel, we indicate /?2- 



23 



o 

00 



5f 



o 

CD 



CD 

0Q 



1 - 



-i— i— i— i— i— i— i— i— i— r 

90-, / \90 K |, 

Zr(n,p) Nb 



- ■ » 



-j— i— i— i— i— |— i— i— r 

(a) Exp 



f 



1 - 



4 
2 




10 



20 



30 



40 



50 





1 1 I 1 1 1 1 I 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 
(b) QRPA - 


. I 


. II I.J , 1 1. 





*f — 1 I — i — i — i — i — i — i — i — i — i — i — r 

10 20 30 



T 



T 



T 



r 

40 



50 



(c) p =0.1 



10 
n I 



20 
I 



I i i i r 

30 40 
— I— i— i— i— i— I— i— i— i— i- 

(d) p = - 0.1 



50 



' 1 I 

10 



~ I I | I I I I | I I I I | I I I I 

20 30 40 50 



E ex (MeV) 



FIG. 6: (Color online) Gamow- Teller strength distributions B(GT + ) on 90 Zr as a function of the 
excitation energy E ex w.r.t. the ground state of 90 Zr. The filled (red) circles in the uppermost 
panels are deduced from the 90 Zr(n,p) reaction at 293 MeV In each panel, we indicate fo. 
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FIG. 7: (Color online) Gamow- Teller strength distributions B(GT ) on 92 Zr as a function of the 
excitation energy E ex w.r.t. the ground state of 92 Zr. Experimental data denoted as filled (red) 



points in the uppermost panels are deduced from the 92 Zr(p,n) reaction at 26 MeV 21]. In each 
panel, we indicate fa- 
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